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ABSTRACT 

We propose a simple model for dark energy useful for comparison with observations. 
It is based on the idea that dark energy and inflation should be caused by the same 
physical process. As motivation, we note that Linde's simple chaotic inflation V = 
T,m 2 (f> 2 produces values of n s = 0.967 and r — 0.13, which are consistent with the 

WMAP la error bars. We therefore propose V = \m\4>\ + \m\4>\ with mi ~ 10~ 5 
and rri2 ^ 10~ 60 , where c = 1 = h and the reduced Planck mass is set to unity. The 
field (f>i drives inflation and has damped by now (</>i,o = 0), while <f>2 is currently rolling 
down its potential to produce dark energy. Using this model, we derive the formula 
8w(z) = w(z) + 1 = 5wo(Ho/H(z)) 2 via the slow-roll approximation. Our numerical 
results from exact and self-consistent solution of the equations of motion for </> 2 and 
the Friedmann equations support this formula, and it should hold for any slow-roll 
dark energy. 

Our potential can be easily realized in N-flation models with many fields, and is 
easily falsifiable by upcoming experiments — for example, if Linde's chaotic inflation is 
ruled out. But if r values consistent with Linde's chaotic inflation are detected then 
one should take this model seriously indeed. 

Key words: cosmological parameters - cosmology: theory, cosmology: dark energy, 
equation of state, inflation, early Universe observations 
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1 INTRODUCTION 

Measuring dark energy is one of the most exciting prob- 
lems in cosmology today. There are a number of expen- 
sive programs underway to measure w, the ratio of the 
pressure to the energy density in dark energy. If dark en- 
ergy is a pure cosmological constant, then w = — 1 for all 
time. Currently, measurements of w are compared to a toy 
model in which w changes linearly with expansion factor a: 
w — wo + w a (l — a). Observational programs are judged 
by a figure of merit which includes their ability to measure 
the quantities wo and w a in this toy model, also known as 
the Chevallier-Polarski-Linder (CPL) parametrization (Al- 
brecht et al. 2006, Chevallier & Polarski 2001, Linder 2003). 

If one parametrizes w(a) in terms of the above toy 
model, the current la limits from the 7-year WMAP data 
combined with BAO+-ffo+SN are wo = —.93 ± .13, and 
w a = —.41 ± .72 (Komatsu et al. 2010). These values are 
consistent at the la level with wo = —1 and w a — 0, which 
would be a pure cosmological constant. 
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In the near future, the Sloan III survey should measure 
the Baryon Acoustic Oscillation (BAO) scale to high ac- 
curacy using LRG's (Luminous Red Galaxies), and should 
lead to a measurement of wo to an accuracy of 3% (SDSS-III 
2008). Using this same dataset, we can use genus topology 
to measure wo independently to an accuracy of 5% (Zunckel, 
Gott, & Lunnan 2010, Park & Kim 2009). Supernova studies 
can achieve similar results (c.f. Riess et al. 1998, Riess et al. 
2007, Albrecht et al. 2006). 

In the longer term, the Euclid and W FIRST (Wide- 
Field Infrared Survey Telescope) satellite missions may be 
able to achieve an accuracy of 1% (Cimatti et al. 2009, 
Blandford et al 2010). Such observations may continue to 
point to wo = — 1 and w a — with higher and higher accu- 
racy, which would be an important result, but would leave 
us still in the dark as to the exact nature of dark energy. 

More exciting would be if a detectable difference Swo 
between wo and —1 is found (i.e. Swo = too + 1). For this 
reason, it is desirable to consider simple models, consistent 
with current data and falsifiable in principle in the near 
future, in which there is a chance that such a detectable 
difference may be found. Such models may be used as a 
guide for interpreting the observations. We therefore propose 
such a simple model, based on the idea that inflation and 
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dark energy come from the same physical mechanism. This 
hypothesis is not present in most of the standard theories 
popular today. 

Today, the most popular theory for dark energy is that 
a string landscape exists with many metastable vacua with 
different values of Vq. In this picture, we are currently sit- 
ting at the bottom of a potential well whose low point has a 
vacuum energy density of Vq. Thus, the accelerated expan- 
sion we see in the Universe today is attributed to our current 
static location at a metastable well in the potential, while in 
contrast, we have evidence that the accelerated expansion we 
see in the early Universe is due to slow-roll inflation, where 
the field is slowly rolling down a potential hill. 

Furthermore, for dark energy today, there is supposed 
to be a complicated potential that is a function of many 
fields fa rather than merely the one field in slow-roll infla- 
tion. The kinetic energy in these fa fields has damped out 
and they are no longer changing. A problematic feature of 
this model is the value of V : 1CT 120 (with c = 1 = 8ttG = ft), 
an extraordinarily small number. The common solution to 
this problem is to propose that there are of order 10 500 vac- 
uum states with values of —1 < Vo < 1. Then one invokes 
anthropic effects to argue that it would be difficult for intel- 
ligent life to evolve unless V ~ 1CT 120 (cf. Weinberg 1987, 
Vilenkin 2003). This is not as satisfying to some physicists 
as an actual prediction of the amount of dark energy we 
observe made directly from the physical model. 

A final problem with this model is the appearance of 
Boltzmann brains (see discussion in Gott 2008 and refer- 
ences therein). Briefly, in the far future, the Universe comes 
to a finite Gibbons and Hawking temperature T = l/(27rro), 
where ro = (3/Vb) 1 ^ 2 , and Boltzmann brains appear. While 
this problem may be manageable depending on what mea- 
sure one uses (cf. deSimone et al. 2010), it is still a problem 
that must be addressed in the popular Vo model. 

This paper will be structured as follows. In section 2, 
we discuss previous models of dark energy and attempts to 
unify dark energy with inflation, as well as observational ev- 
idence for inflation. In section 3, we discuss Linde's chaotic 
inflation. In section 4, we propose our double inflation poten- 
tial based on Linde's model. In sections 5 and 6, we show how 
our potential might be realized as an effective N-flationary 
potential arising from many axion fields, and explain why 
this may be regarded as preferable. In section 7, we esti- 
mate the probability of different deviations from w = — 1. 
In section 8, we describe the numerical approach we used 
to obtain exact values of Swo and w(z), and in section 9 we 
present the results we obtain. In section 10, we conclude. 

We use c — 1 — h = 8-kG throughout, where G is New- 
ton's constant. In these units, the reduced Planck mass m p i 
= (hc/8TvG) 1/2 = 2.44 x 10 18 GeV is equal to unity. These 
units follow Linde (2002). 



2 PREVIOUS DARK ENERGY AND 
UNIFICATION MODELS 

2.1 Dark Energy 

While a pure cosmological constant w = — 1 is the simplest 
model that explains current observations (see Li, Li & Zhang 
2010), there are numerous other dark energy models, which 



we briefly discuss here to place our own approach in context. 
Broadly, there are scalar-field models such as quintessence, 
k-essence, a tachyon field, phantom dark energy, dilatonic 
dark energy, and a Chaplygin gas. Alternatively, there are 
also proposals that the Universe's accelerating expansion 
may stem from modified gravity, effects of physics above 
the Planck energy, or backreactions of cosmological pertur- 
bations. For these models, we direct the reader to the review 
of Copeland et al. (2006); here we will focus on summariz- 
ing the scalar-field models, since our model falls into this 
category. 

Quintessence is characterized by a scalar field that rolls 
down its potential to produce the accelerated expansion we 
observe. Slow-roll solutions are possible if certain conditions 
are satisfied (see Copeland et al. 2006), and both exponen- 
tial and power-law type potentials are popular. One possible 
difficulty with quintessence models is that if the field cou- 
ples to ordinary matter this could lead to time dependence 
of the constants of nature. In contrast to quintessence, for 
k-essence, the expansion stems from the kinetic energy in 
the field (Armendariz-Picon et al. 2001, 2000, 1999). These 
models have non-canonical kinetic energy terms. 

In addition, there are tachyonic fields, which have been 
a candidate for both inflation and dark energy. For further 
discussion, we refer the reader to Copeland et al. 2006, Pad- 
manabhan 2002, Bagla et al. 2003, and Abramo & Finelli 
2003; the most notable feature of the tachyon models is that 
the equation of state varies between and —1 regardless of 
how steep the tachyon potential is; for accelerated expan- 
sion, however, one must have a shallow potential compared 
to V(fa) oc <t>- 2 . 

The models enumerated above all correspond to w > 
— 1; however, w < — 1 has not been observationally disal- 
lowed, and this region of parameter-space is generally pro- 
duced by so-called "phantom" or "ghost" dark energy. These 
are scalar fields with negative kinetic energy, and suffer from 
a number of difficulties, most notably quantum instabilities 
(see Copeland et al. 2006 and references therein). Dilatonic 
models are attempts to stabilize phantom dark energy (Cline 
et al. 2004, Gasperini et al. 2002). 

Finally, there is the Chaplygin gas, with equation of 
state p oc —p, which at early times behaves as a pres- 
sureless dust and at late times as a cosmological constant 
(Copeland et al. 2006). It is thought that the Chaplygin gas 
would result in a strong Integrated Sachs- Wolfe effect, and 
hence loss of power in CMB anisotropics; though this can 
be avoided, the constraints are rather tight (Kamenshchik 
et al. 2001). 

2.2 Inflation 

Since our model is based on the idea that inflation and dark 
energy are the results of the same physical process, we briefly 
review the evidence for inflation here to contextualize our 
model. Proposed by Guth in 1981, inflation, the idea that the 
Universe underwent a period of exponential, super-luminal 
expansion in the first 10 -35 seconds after the Big Bang, 
explains many "puzzles." For instance, the isotropy of the 
CMB down to one part in 100, 000, even comparing regions 
that are not causally connected in the standard Big Bang 
model, is a natural feature of inflation (see Komatsu et al. 
2010, Guth 2007). 
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Further, the Universe is flat down to a few percent (see 
Tegmark et al. 2004); Dicke and Peebles (1979) point out 
that to realize this flatness today in the standard Big Bang 
model, fltot would have had to be unity to 59 decimal places 
at the Planck time, an extraordinary fine-tuning. But in in- 
flation, the entire observable Universe today came from such 
a small initial region in comparison to the radius of curva- 
ture of the early Universe that flatness arises generically and 
requires no fine-tuning. 

Another piece of evidence for inflation is the absence 
of magnetic monopoles. All Grand Unified Theories pre- 
dict them. Indeed, Preskill (1979) estimates they would be 
dominant in a non-inflationary cosmology. However, none 
are observed. Inflation provides a natural mechanism to ex- 
plain why, as the exponential expansion moves the nearest 
monopoles beyond our observable Universe. Finally, infla- 
tion predicts a density spectrum of almost scale-invariant, 
adiabatic Gaussian fluctuations, a spectrum that fits the ob- 
served temperature fluctuations in the CMB as a function 
of angular scale quite well (Guth & Kaiser 2005, Guth 2007, 
Komatsu et al. 2010). 

Thus it appears we experienced an epoch of accelerated 
expansion in the early Universe (inflation), and we observe 
that we are in another epoch of accelerated expansion today 
due to dark energy. With all of this in mind, it seems the 
simplest model of dark energy would be one in which it 
and the inflation we encounter in the early Universe are 
essentially identical. 



2.3 Unification Models 

Unification of dark energy and inflation has been pro- 
posed before. For instance, Peebles and Vilenkin (1998) 
used a scalar field potential V oc (f> 4 + m 4 , (f> < 0, V oc 

8 

tpi+mi i 4> > 0, which is supposed to produce inflation like 
a self-interacting field for <f> < and then dark energy as it 
slow-rolls for cf> > 0. Cardenas (2006) proposed a tachyonic 
field with an exponential potential V oc exp[— (j>/(f>o], where 
there is inflation followed by reheating and finally an epoch 
of dark energy inflation as 4> continues to roll down the hill of 
its potential. Finally, Liddle & Urena-Lopez (2006) explore 
two scenarios: unification of dark matter, dark energy, and 
inflation into the same scalar field, and unification of dark 
energy with dark matter in one scalar field while inflation is 
provided by another. 

Double inflation scenarios in which the early Universe 
undergoes two separate epochs of inflation have also been 
proposed (these scenarios do not attempt to explain dark 
energy as coming from an inflaton-like field). For instance, 
Poletti (1989) considers the quantum cosmology of a mini- 
superspace model, and focuses on the idea that the Wheeler- 
DeWitt equation may be invariant under rescalings of the 
metric. We note this work simply to give credit to Poletti 
for the idea of a double-inflationary scenario (he also derives 
the relations showing how the two fields' evolutions relates 
to their masses) ; his focus on the Wheeler-DeWitt equation 
is very different from our own on dark energy. 

Turner et al. (1987) suggest double inflation in the early 
Universe as an exlanation for a bimodal density fluctuation 
spectrum (fluctuations on small and large scales) consistent 
with the observations available in the late 1980's on galaxy 



structure. Silk & Turner (1987) make a similar suggestion, 
arguing that the decoupling of large and small scale struc- 
ture such a scenario would produce could explain the (circa 
1986) excess of observed power on large scales. Muller and 
Schmidt (1988) propose that double inflation might arise 
from fourth-order gravity coupled to a scalar field, and sug- 
gest that the second epoch of inflation would yield dark mat- 
ter. 

Finally, Polarski & Starobinsky (1992) calculate the 
spectrum of perturbations produced by double inflation, and 
are again motivated by the apparent mismatch between the- 
ory and observation in prediction of large-scale structure. 
More recently, Yamaguchi (2001) uses double inflation to 
explain large-scale structure (first epoch) and small-scale 
structure (second-epoch), suggesting it might naturally arise 
from supergravity; Yamaguchi also does not suggest double 
inflation as an explanation for dark energy. 



3 CHAOTIC INFLATION 

Let us begin with Linde's (1983) chaotic inflation, arguably 
the simplest model of inflation ever proposed. Linde's po- 
tential was of the form 



V(cf>) = \m 2 <f 



(1) 



in Planck units where c = 1 = h = 8nG, m p i = 
(Hc/8-kG) 1 ^ 2 = 1 is the reduced Planck mass, and G is New- 
ton's constant. 

This is a simple massive scalar field. The full equation 
of motion for the field is 



• + 3Hcj> - 



dV 

~d(j> 



(2) 



where H is the Hubble constant and 3H(fi is a frictional 
term due to the expansion of the Universe (see Copeland 
et al. 2006, Linde 2002, Lyth & Liddle 2000). In general, in 
inflationary models where the Universe is effectively flat, 



W 2 1 
H =-p 



(3) 



(with SttG — 1). Taking p to be due primarily to the 
inflationary potential V(4>), we have 



H = 



1/2 



(4) 



The approximation follows by noting that in slow-roll infla- 
tion, the condition i\? << V(4>) is satisfied (see Copeland et 
al. 2006, Linde 2002). Thus p = \tf + V{4>) w V{4>). 

We now make the further approximation that cf> « <fi, 
which is satisfied because the frictional term 3H<f> means 
that the field quickly reaches a nearly constant "terminal" 
velocity. We thus have from equation (2) that, in the early 
Universe during the inflationary epoch, 



m 2 <f> 

Since H — dlna/dt, we have 



d In a = H 



dt_ 

d4 c 



(5) 



(6) 
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where the last equality follows by using equation (4) for 
H and equation (5) for dt/d(j) = 4> _1 . 

If N is the number of e-folds of inflation then 



N = In 



("final 



^initial 



(7) 



(^initial 

Inflation continues in slow-roll until (j> 2 finai ~ lj when 
the kinetic energy | if 2 becomes comparable with the poten- 



tial energy 



,2/2 



At this point, the exponential expansion 



ends and the kinetic energy in the field is dumped into the 
thermal energy of particles, ushering in the hot big bang 
epoch — a radiation dominated, thermal- energy filled Uni- 
verse. The (f> field is damped and settles at <j> — during this 
process, so that V(4>) = and the vacuum energy density in 
the massive scalar field becomes zero. (Or the tiny value of 
v = icn 120 if one adds a tiny constant to the formula such 
that V{4>) = Vo + ^m 2 (f) 2 to account for dark energy. But 
we will not be adding this Vb term.) 

Fluctuations re-entering the causal horizon today left 
the causal horizon approximately N = 60 e-folds prior to the 
end of inflation, when according to equation (7) the value of 
<j) 2 = 4iV + 1 = 241. In this case, the value of the power law 
primordial tilt evaluated at ko = .002 Mpc -1 should be 



1+2 



V 



1 A» 



.967 (predicted) (8) 



(cf. Easther and McAllister 2006) and the value of r, 
the ratio of tensor to scalar modes, should be 



^2 ~ n 



.13 (predicted) 



(9) 



(cf. Kim and Liddle 2006) . Fitting the amplitude of the 
observed fluctuations requires 

m = 7.8xl0~ 6 (10) 

(cf. Kim and Liddle 2006). Remarkably, the predicted 
values of n 3 and r are consistent with the observed values 
from WMAP+BAO+ffo (Komatsu et al. 2010): 



n s = .968 ± .012 (observed) 

and the 95% confidence level constraint 
r < .24 (observed). 



(11) 



(12) 



The agreement between the predicted and observed val- 
ues of n 3 is especially impressive considering that potentials 
of the form V((f>) = (A/4)0 4 have been ruled out by pre- 
dicting unacceptable values of n s — .95 and r = .26 (cf. 
Komatsu et al. 2010). Given this, the search for the ten- 
sor modes (r > 0) is on — for instance, the Planck satellite 
hopes to improve the measurement of r. Polarization studies 
in the future should if successful offer a smoking-gun proof 
that the tensor modes are there. Such modes are not pre- 
dicted by the Ekpyrotic/Cyclic scenario and, if found, they 
would offer a convincing proof of inflation (cf. Linde 2002). 

It is remarkable that a model as simple as Linde's 
massive-scalar-field chaotic inflation is currently consistent 
with the observational data. Linde's model predicts in par- 
ticular a value of n s noticeably less than 1 — a hallmark of 
slow-roll inflation that is indeed observed. 

While more complicated potentials can produce lower 
values of r (cf. Kallosh and Linde 2010), this comes at the 
expense of adding more free parameters. The Linde theory, 



in contrast, is simple enough that it offers the possibility of 
being confirmed in a dramatic way if the observed value of 
r is .13. If that occurs, we will no doubt conclude that the 
inflation seen in the early Universe is due to a massive scalar 
field. In that case, we argue here that we should expect a 
similar origin for dark energy as well. 



4 DOUBLE INFLATION 

4.1 Swq in terms of the dark energy scalar field 
today 

We propose the following simple double-inflation potential 
for inflation and dark energy: 

T/ 1 2,2 . 1 2 ,2 

V = -m^i + -m 2 <j> 2 



(13) 



with mi ~ 10~ 5 and m 2 ^ 10~ 60 . Double inflation 
with a potential of this form was introduced (typically with 
mi ~ (20 ± 5)m,2) to explain inflation alone (Polarski & 
Starobinsky 1992, see also Silk & Turner 1987). We will be 
using it with widely different mass scales to explain inflation 
and dark energy. The equation of motion for the inflaton 
field is 



k + 3Hcf>i + m 2 (f>i = 0, 

and for the dark energy field is 



0. 



(14) 



(15) 



As noted earlier, we are working in Planck units where 
c — 1 — h — 8ttG, m p i = (fic/87rG) 1//2 = 1 is the reduced 
Planck mass, and G is Newton's constant. 

In the inflationary epoch when the Universe is dom- 
inated by the vacuum energy density provided by V, we 
have 

3H 2 = V = ^m 2 2 + ^m 2 2 4> 2 2 . (16) 

The slow-roll approximation is valid for both fields and 
the evolution of the two fields is given by 



m- 



02 , initial 



mi 



In 



L(i) 



(17) 



(cf. Easther & Mcallister 2006). Since m2 << mi, 
<p 2 (t) w constant w cf>2,initial even though <t>i(t) / $1, initial 
evolves considerably during inflation. There have been 60 e- 
folds of inflation since the perturbations now re-entering the 
horizon left the causal horizon, but there could have been 
more e-folds of inflation before that, so we expect iV > 60 
and thus \[(t>l, in itiai-<t>\jinai] = N > 60. Since <j)\j irla i = 1 
marks the end of inflation, 4>\ tin m a i > (4 x 60 + 1) = 241. 
As inflation ends, the kinetic energy in the <j>i field is con- 
verted into thermal particles, the motion in <f>i damps, and 
<f) 1 comes to rest at a value of cf>\ = 0. The dark energy seen 
today thus derives from the 02 field, and V — \m%4>% today 
because 0i.o = 0. 

Since we observe a dark energy acceleration today con- 
sistent with w rs — 1 we expect slow-roll inflation to apply 
|j<-i| << 3H\<fr\, <j} 2 « V(<£)j. Since dark energy is domi- 
nant today, we use our results from Section [3] Since today 
Pde.q = Q.DER0, po = (1/Qde)pde,o- Making these changes 
to equation (4), we have 
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H 



-m<2,4>2fl, 02,0 



3H 



2Qr 



-m 2 .(18) 



Now, w, defined to be the ratio of pressure to energy density 
in the dark energy, is given by 



P 
P 



(19) 



where the equality is specific to scalar quintessence 
fields (see Copeland et al. 2006) , and we recall that because 
the inflationary field has damped long ago, V = \m\<j>\ in 
the present epoch. We therefore have from equations ( 18 1 



and (191 that 



W : 



*2,0 



The difference between wo and 
be 8wq, is thus 



Swq = UlQ + 1 



(20) 

-1, which we define to 
(21) 



Current limits from WMAP suggest Swq = .07 ± .13, 
where 5wo is simply Sw today. Importantly, since in the 
future we expect 02 to roll down to zero, leaving a vacuum 
energy density of zero (with no Vo term), the Boltzmann 
brain problem disappears. See Figure [2] for a comparison of 
this formula to our numerical results from a full and self- 
consistent solution of the equation of motion for 02 and the 
Friedmann equation as described in Section[8] Equation (211 
is only valid evaluated today, because today we are dark- 
energy dominated, which is the assumption under which this 
formula was derived. This formula will not be valid as a 
function of time in the past, because dark energy has only 
recently become dominant. 



4.2 <5w as a function of redshift 

In what follows, we will derive a formula that gives Sw as a 
function of time. Let us be very clear here that we are no 
longer assuming dark-energy domination as we did in |4.1| 

The formula we present below will give 5w(z) as a func- 
tion of H(z) and so is a simple template observational pro- 
grams can use to look for slow-roll dark energy. 

We have from equation ( |19[ ) and the definition Sw — 
w + 1 that 



Sw = 



(1/2)0| + V 



(22) 



From this, we can see that small 8w implies that 2 << V, 
and so we approximate that 



5w 



0_f 
V 



(23) 



Since 02 is small, (j>2 does not vary much with time and so V 
is approximately constant in time. Thus 5w tx 2 . S ince in 



slow-roll 02 is also small, we have from equation ( 15 \ that 



37702 



dV 



= -m 2 2 - 



As we have already noted, 
Hence 3H<j) 2 ~ constant, so 



(24) 

62 is approximately constant. 



t>2 OC 



Sw OC 



H 

Since Sw oc 02, equation (25) implies that 
1 



H 2 ' 



Normalizing appropriately, we obtain 



8w(z) ~ 5wo 



Hp 

H{z) 



(25) 



(26) 



(27) 



In summary, then, for small Sw, slow-roll applies, 2 is 
small, and Sw oc 2 . When the field reaches terminal veloc- 
ity, the acceleration 2 is nearly zero, so the Hubble friction 
term 3/f 2 is balanced by the slope of the potential the 
field is rolling down, —dV/d(/>2- This is analogous to a ball 
rolling down a hill: it will reach terminal velocity when en- 
ergy dissipation from friction cancels out the energy it gains 
by moving to lower values of its potential. Finally, if 02 is 
roughly constant in time, dV/d(f>2 will be as well, meaning 
02 oc H~ , which leads to equation (26). This derivation 
could easily be duplicated for other dark energy potentials 
as long as Sw is small. 



We term this formula (equation ( 27 1 ) our "slow-roll for- 
mula" in subsequent discussion. See Figure [3] for a compari- 
son of this formula to our numerical results from a full and 
self-consistent solution of the equation of motion for 2 and 
the Friedmann equation as described in Section [8] 



5 N-FLATION 

A potential of the form in equation (13) governed by equa- 
tions of motion (14) and (15) can be produced easily by 
models of N-flation. A possible criticism of the original Linde 
chaotic inflation is that it requires > 1 in Planck units. 
Linde argued that this was acceptable as long as V < 1. 
But it was felt that it would be difficult to produce values of 
> 1 in string theory models. Thus, N-flation (Dimopoulos 
et al. 2008) has been proposed (cf. also Easther and McAllis- 
ter 2006), motivated by the fact that supersymmetric string 
theories allow of order 10 5 axion fields. 

For such axion fields V(ipi) = /i 4 [l — cos(-0i//i)], and 
for ipi « 1 (i.e. tpi significantly below the Planck mass — 
which we would like) we note that the potential is of the 
Linde quadratic form with effective mass m — u 2 / fi. As we 
have already pointed out, string theory allows the number 
of such fields to be large (Dimopoulos et al. 2005). Hence in 
this work we will adopt N = 10 4 . 

In this model there are JV fields tpi (where i = 1, . . . , N) 
with approximately equal masses « m. Then the poten- 
tial is V = £ \ra}^} » 1/(0) = |m 2 2 where 2 = ]T> 2 . 
This is so-called "assisted inflation." 

Since each rrii « m, all the ipiS evolve together via 



In 



yji,initial 



In 



^Pj, initial 



hi 



1pj,initial 



(28) 



if nrii w nij for all i and j. This will be true if the mass 
spectrum of the fields is strongly peaked and densely packed: 



i 2 = m 2 exp [(J, — l)/er] 



(29) 



where a > 280 for N > 600 (Kim and Liddle 2006). 
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If the fields are strictly non-interacting the masses could in 
fact be exactly equal. 

We want double N-flation with a potential 

V = ^2 imfVi.i + ^aV'JM = ^"^101 + ^ m 202 (30) 

where there are N = 10 4 ip 2 fields each of mass m2 10 
and N = 10 4 ipi fields each of mass mi ~ 10 -5 (in keeping 
with our hypothesis that inflation and dark energy should 
arise from the same process) and by definition (j>2 ~ V , 2,i 
and analogously for cf>\. Initially we need (pi initial > 241 
to produce enough inflation (> 60 e-folds) to explain our 
Universe. With 10 4 fields, that just means that each ip 2 ^ > 
.0241 and so all the ipij's can be sub-Planckian (< 1). This 
is good. 

Are such low values of mj w 10 -60 plausible from 
the point of view of string theory? Interestingly, Kaloper 
and Sorbo (2006) have independently proposed just such 
an N field quiNtessence model for dark energy using ul- 
tralight pNGB [pseudo-Nambu Goldstone bosons] (axions) 
from string theory. They argue for potentials of the form 
V(ipi) = fi 4 [l — cos[ipi/fi)]. Svrcek (2006) has also argued 
for multiple ultra-light axion fields with potentials of this 
form to explain dark energy. 

We note in each case that for sub-Planckian ipi, initial's, 
the potential is of the desired Linde quadratic form with 
effective mass m = /i 2 /fi- Svrcek notes that pseudoscalar 
axion fields have a shift symmetry and if this symmetry were 
exact it would set the potential to zero and the axions would 
be massless. In string theory the shift symmetry is broken 
only by nonperturbative effects. In string theory axions thus 
receive potential only from nonperturbative instanton effects 
which are exponentially suppressed by the instanton action. 
Hence, if the instantons have large actions they can give rise 
to a potential many orders of magnitude below the Planck 
scale. 

Svrcek argues that /i 4 = M 4 exp[— Si ns t], where M ~ 1 
and Sinst ~ 280, can create a vacuum energy density today 
comparable with what we observe for dark energy. Svrcek 
adds a Vb term as well, which we eliminate as unnecessary. 
We argue that if the axion fields are able to explain the 
amount of dark energy we observe today the Vb term can 
be eliminated. Both Kaloper and Sorbo (2006) and Svrcek 
(2006) are explicitly creating quintessence models for dark 
energy. Both also note that single field models with sub- 
Planckian field values are unacceptable for quintessence and 
favor models with N = 10 4-5 fields. 

We are proposing to combine these quintessence mod- 
els that use N fields with the N-flation models to explain 
dark energy and inflation. Independently Svrcek (2006) also 
speculates "Hence, it could be that some of the string the- 
ory axions have driven inflation while others are currently 
responsible for [a] cosmological constant." We take the point 
of view here that there are N = 10 4 equal-mass ultra-light 
fields ip2,i that create a slow-roll dark energy and N = 10 4 
equal-mass heavy fields ipi^ that create slow-roll inflation in 
the early Universe. 

If there are in addition singleton fields with intermedi- 
ate masses, with sub-Planckian ipi's also, these would not 
have inflated but rather would have rolled down, as Svrcek 
notes. Some of these axion fields could have rolled down 
and created dark matter. They do not cause inflation be- 



cause ipi < 1, so when the other thermal particles red- 
shift so that the ipi vacuum field energy becomes dominant, 
3H 2 = \m\ip\. Thus H is not large enough to cause the low 
velocity (ipi) required for slow-roll inflation. 

However, if there are many fields of essentially the same 
mass, 3H 2 = ^mfVl ~ § m *V 2 > where (p 2 > 1. In this 
case, H is much higher (by a factor of %/N) causing each ipi 
to be lower and slow-roll inflation to occur. Thus, inflation 
only occurs when many fields congregate at the same mass 
scale. 



6 WHY N-FLATION IS SUPERIOR 

We expect all of the ip2,iS and ipi^'s to be sub-Planckian. 
If N = 10 4 , this means that (pi initial < 10 4 , which means 
(using equation 7) that there can be at most 2500 e-folds 
of inflation in our Universe. Hence our Universe today is 
less than exp[2440] times larger than the visible horizon. 
In Linde's original formulation of chaotic inflation, random 
quantum fluctuations allowed Universes to give birth to Uni- 
verses with various values of <pi,imtiai- The ones with larger 
values of <p>\. initial grew faster until most of the volume of 
the multiverse was in the fastest expanding states, with 
V = 1 = \m\(p\. That would mean (pi initial ~ 10 10 and 
the Universe today would be exp[10 10 ] times larger than the 
part we can see. 

For our model of dark energy, a simple Linde-type 
chaotic double-inflation picture would eventually lead to 
most of the volume of the multiverse being in the fastest 
expanding states, given by the ellipse 

V = 1 = -mUi + \ml<t>l 

= & + & (31) 

2 x(10 5 ) 2 2 x(10 60 ) 2 ' 

This ellipse is very elongated in the <p2 direction, and at 
random points on it the <p2 field contributes just as much to 
the potential and to the inflation as the (pi field. Thus start- 
ing values <p2,initiai ~ (mi /m,2)<pi, initial would be expected, 
and since (p2 is slower to roll down than (pi, we would not 
get the sub-dominant dark energy we require. 

If we use N-flation to realize the potential in equation 
(13) via equation ( |30[ ) this problem does not occur. Since all 
of the fields are sub-Planckian, the fastest inflating regions 
are characterized by starting values of cp2 ~ (pi that are 
bounded above by N, and since m,2 « mi, the contribution 
of the (p2 field to the potential is sub-dominant (< \m\N 
for the (p2 field versus < ^m 2 N for the (pi field). 



7 BOUND ON PROBABILITIES FOR 
DEVIATION OF w FROM 1 

In this picture we might expect the initial values for (p 2 and 
(pi to be comparable. What is the smallest <pii n iu a i could 
be? It must be at least 241 to explain the at least 60 e- 
folds of inflation we see within the visible Universe. By the 
above argument, we might expect (p2,initiai to be similar. 
Because of equation (17), (p2 does not evolve much during 
the period of inflation. Its main chance to roll down is in 
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above. However, the number of galaxies (and observers like 
ourselves) produced in our Universe is proportional to 



exp[3iV] = exp 



^1 , initial 



(33) 



So it might not be surprising for us to observe 
< <f>i,initial > less than < 4>\,imtiai >, but the details depend 
on questions of measure which are unsettled. 

Suffice it to say, one must find a measure that makes 
what we observe, namely 0| initia i > f 2 (since Swo = 7% ± 
f3% from WMAP) and 4>\ >iniUa i > 241 (to produce at least 
60 e-folds of inflation) , not particularly unlikely. For the time 
being, we suggest taking an empirical approach and asking 
what future observations can tell us about 02 initial. 



Figure 1. Plot of our upper bound on the probability that 5wo is 
greater than a given value (equation j32[ l). This is an upper bound 
because in fact (pi initial could be much greater than 16 (chosen to 
provide approximately the minimal 60 e-folds of inflation), leading 
to larger values of 02 and lower values of (5u>o. However, the plot is 
moderately encouraging: there still could in principle be a nearly 
20% chance that Swo could be as large as 3%, which would likely 
be detectable. 



the current epoch when H is low. But there have not been 
many e-folds of inflation during the current epoch, so we 
might expect 0|,o to be only a little less than its minimal 
initial value of about 241. That would give a value today of 
Sw ~ 4/ ((30! >o /O D s) + 2) = .4% via equation (21). 

Since we expect the 02,i's and the ^>i,j's of equation (30) 
to be uncorrelated in our n-flationary picture, we expect by 
the central limit theorem that 02, initio; and (/>i,initiai are 
Gaussian. If their magnitudes are comparable, in keeping 
with our hypothesis that the physical processes for inflation 
and dark energy are identical, we might a priori expect on 
average to find < 4>l Anitial > = < <j>\, in itiai > ■ 

Consider the probability function 

P($2,initiali4>\,initial)<i<i>-l,initial<l<i>\,initial- If both Vari- 
ables follow a Gaussian distribution, P has circular 
probability contours in the (0 2 , initial, </>i,initiai) plane, and 
the probability of finding \<j)2,initiai\ < X|0 Mniti( a| for 
X < 1 is simply P = (2/7r)arctanX. If <j>\ >initial > 241, the 
probability of observing a value of Swo > Sw' is 



P(8wq > 5w' ) < — arctan 

TV 




IflDE 


4 


- 2 




Jw' 





(32) 



using equation (21) to relate 02, o — 4>2,initiai to Swo. See 
Fig. 1. The probability of observing 8w > 3% is P < 20%. 
It could be considerably smaller if there are significantly 
more than the minimum 60 e-folds of inflation in the early 
Universe. 

With this in mind, it may be reasonable to expect a 
small value of Swq. On the assumption that the same phys- 
ical processes led to both inflation and dark energy, if there 
were at least 60 e-folds of inflation in the early Universe, 
there should be at least 60 e-folds of inflation due to dark 
energy ahead in the future, as indicated in the calculation 



8 NUMERICAL METHOD 

To integrate the equation of motion for 02 (15) numerically, 
we must also solve the Friedmann equation for H as a func- 
tion of time. Since we know from observational constraints 
that Sw is small, we start by assuming that w = — 1 and 
solve the Friedmann equation to determine H as a function 
of time. Then using this H we can solve the full equation of 
motion for 02 as a function of time (though we will instead 
solve the re-scaled equation (43) as described later in this 
section). Knowing 02 gives us w as a function of time, and 
we can insert this back into the Friedmann equation to give 
us an improved H as a function of time. 

We then iterate. This approach converges rapidly on 
a self-consistent solution where we have w, H, and 02 as 
functions of time and both the Friedmann equations and 
the equation of motion for the field 02 have been solved 
self-consistently and exactly. 

For starting conditions at t = we assume d(f>2/dt=0 
since H(t) tends to infinity as t tends to zero. We discuss 
the starting condition on 02 itself later. Since in this sec- 
tion we will consider only the 02 field with mass m2, for 
greater legibility we suppress the subscripts on 02 and rri2 
in what follows. We recall equation (15) (with subscripts 
suppressed): 



3.ff0 + m 2 = O, 



(34) 



where H — (l/a)(da/dt) is the Hubble constant. Defin- 
ing a new variable for time r = t/tno = tHo and transform- 
ing and 0, we have = 4>"Hq and = 0'i/o, where prime 
denotes a derivative with respect to the new time variable 

T. 



We thus have from equation (341 that 
,H(t) 



+ 3- 



H 



+ 



H- 



= 0. 



(35) 



and its derivatives are functions of r, so we desire that 
H should be as well. We convert the Friedmann equation 
from time variable t to time variable r; so doing will yield 
the Friedmann equation governing H(t) = H(t)/Ho- 

It is 



H 1 



+ Q,de exp 



3^ T °/7(r') [(l + w(r')]dr' ), 



(36) 
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where Q r = p r /pc, and analogously for fi m and SIde- 
p c = 3Hq/(8wG) is the critical density; subscript r denotes 
radiation, m matter, and DE dark energy, to is the value of 
r today (when a — 1). Note that for w = — 1, this reduces to 
the usual Friedmann equation used in ACDM models. We 
thus find the equation for 



+ 3H(t) 



i 2 

+ 



= 0, 



(37) 



where we have defined m 2 = m 2 /Hq. 

In principle, we can now solve the Friedmann equa- 
tion numerically for H(r) and using it obtain numeri- 
cally. However, there is a constraint on 02, o that we must 
satisfy as we solve numerically. The Universe is flat with 
Q.de + ~ 1. Since dark matter is present, the value of 
Hq is larger than it would be if only dark energy were present 
by a factor of fi^; WMAP-7 gives Qde = -73 (Komatsu 
et al. 2010). Thus at present 



3# = 



2Q. L 



-m 2 (j)lX, 



(38) 



where X = (0o + m 2 0o) /«i 2 0o- 

Using the definitions of w and Sw, we find that 

2 _ 3Qde(2 — Sw ) 

l * J.2 

02,0 



(39) 



Hence we are not free to choose both m* and (j>2,initiai 
independently because 02, initial will determine 02, o, which 



must be consistent with m* such that equation ( 39 1 is satis- 



fied. We therefore require a method of solving the equation 
of motion where we can set 02, initial (and hence 02, o) af- 
ter we already have a solution. This motivates us to observe 
that the equation can be dynamically rescaled by writing 
0'/0 = C( r )- ^ is evident that is always non-zero. Writing 
0' = C0, we find that 0" = (C 2 + C') 4>- Substituting these 
relations into equation (37), we obtain 



C 2 + C' + 3H(t)( + ml = 0. 



(40) 



We choose m» and numerically solve this equation be- 
ginning at r = 0, where we have £ = because of our earlier 
comment on 02,»n««ot- This is done iteratively in conjunc- 
tion with the solution of equation (36) as in the overview 
at the beginning of this section, and in greater detail below; 
the final iteration gives us Swq. We can then determine 02, o 



by inverting equation ( 39 \ 



We have already given an overview of our iterative pro- 
cedure; here we give details. As noted before, we begin by 
finding H as a function of r for w = — 1, which we do us- 
ing the Friedmann equation (36). For w = — 1, this does 
not require knowing to. However, solving for to (= t such 
that a = 1) in this first iteration does provide to for the 
next round of iteration, where we will use w = w(r) in 
equation (36). We use H as provided by (36) in the dy- 
namically rescaled equation (40) for £ = 0'/0, which yields 
w — (£ 2 — m 2 )/(£ 2 + m 2 ) via equation (19). We insert this 
w(t) back into (36) and find a new H, which goes back into 
(40), again yielding w. We continue this process until the 
difference between the n th and n + 1 th step integrated over 
all time is negligible compared to Swq. Then we have an ex- 
act and self-consistent solution to the equation of motion for 
the 02 field and the Friedmann equation. 



(5w today vs. 02 today; 
solid blue is numerical, 
red dashed is approximate 
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Figure 2. Plot comparing the "back-of-the-envelope" approxi- 
mate formula Swo — 4/( j-p — cj>\ q + 2) with the numerical results 
obtained as outlined in Section[8] The plot shows that, while the 
back-of-the-envelope calculation captures the order-of-magnitude 
of Swq and the rough shape of the curve, it should not be trusted 
in lieu of the numerical results. Since it was derived on the as- 
sumption that Swq << 1, it makes sense that the approximate 
and numerical results are closer in this region of the plot (i.e. for 
larger 2 , o) • 



9 NUMERICAL RESULTS 

In this section, we discuss our numerical results. Figure [2] 
compares our approximate formula for Swo as a function 
of the value of the dark energy scalar field today, (equation 
(21 1), to our numerical results from exact and self-consistent 
solution of the Friedmann and scalar field equations, de- 
scribed in Section 8. The curves have similar shapes. How- 
ever, the approximate formula is not a good substitute for 
the full numerical results in the range 4 < 02, o < 12. For 
02,o > 12, the match between analytical and numerical is 
somewhat better. Since our approximate formula was de- 
rived in the limit that Swo « 1, it is not surprising that it 
agrees better with the full numerical results the smaller <5 TO ,o 
becomes. 

Figure [3] compares our "slow-roll formula" (equation 
(27l) with our full numerical results from exact and self- 
consistent solution of the Friedmann and scalar field equa- 
tions. Note that the function H used in these plots comes 
in each case from the self-consistent solution of the Fried- 
mann equation as described in section 8. The "slow-roll" 
formula agrees well with our numerical results and best for 
Swo < 1%. Analogously to our discussion of Figure [2] this is 
not surprising because the "slow roll formula" was derived 
in the same limit that Swq << 1. Hence the smaller 8wq is, 
the better we would expect the agreement to be. 
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<Sw as a function of redshift 
for 02,0 = 16, 10, 6, and 4 
with analytic formula (dashed) 
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Figure 3. Here, we show our numerical results for Sw as a 
function of redshift and also the "slow-roll formula" (dashed) 

Sw = 8u>o ( 1 (see Section 4 for derivation). The values of 
<f>2,0 listed on the plot go 16, 10, 6, and 4 from bottom to top for 
the solid curve/dashed curve pairs. For smaller <5u>o the assump- 
tions underlying our "slow-roll formula" are better satisfied, and 
one can see it agrees more closely with the numerical results in 
this regime. 



w as a function of redshift for 02, o = 6 
and for toy model (upper curve) 

-0.975 , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — r- 



-0.980 



-0.995 




-1.000 - , , , , 
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Figure 4. Here, we compare our model's predictions for w as 
a function of redshift to those of the popular toy model w = 
wq +M>a(l — a), where wq and w a are constants and a is the scale 
factor appearing in the Friedmann equation. The toy model, fit at 
both ends, is the shallower, upper curve. The comparison shows 
that it will be harder to observe deviations from w = — 1 in the 
past if our model is correct than if the toy model is correct. 



Figure [4] compares a representative result from our 
model (with 8wo = 2%) to the Chevallier-Polarski-Linder 
parametrization, w = Wo + W a (l — a). We chose the param- 
eters wo and w a so that the toy model agrees with ours at 
a = 1 and at a = (today and at the Big Bang). This 
leads to wo = —1 + Swo and w a = — Swo- Since the CPL 
parametrization is the shallower curve in the plot, it is clear 
that if our model is correct, it will be harder to observe w 
different from —1 in the past than if the CPL parametriza- 
tion is right. A comparison of Figures [3] and [4] shows that 
our slow-roll formula provides a much better fit for our nu- 
merical results than the CPL parametrization, as it has a 
physical motivation from the slow-roll approximation. 

Figure[3]tells a story consistent with Figure|2] the larger 
02, o is, the smaller is Swo. The main take-away is that 5w 
is a quite steep function of redshift; as the figure shows, 
Sw is nearly indistinguishable from —1 for z > 3 for all of 
the curves (different curves correspond to different values of 
02,o-) This means that if our model is correct, it will be very 
challenging to observe deviations of w from —1 except in 
the recent past. If our model is correct, it will also be more 
challenging to observe deviations from w = — 1 than if we 
assumed that 8w was constant and non-zero. 

Note that the initial values of 02 to which the curves in 
Figure [3] correspond are roughly the same as the final values 
we have quoted in the Figure: (02, initial — 02, o)- This is be- 
cause the Hubble friction ensures that 02 does not roll down 
much (see equation (15)). We chose 02. o = 16 on the as- 
sumption that 02,0 — 4>2,initial — <f>i,initial an( l 01, initial —16 

because that is the minimum value that can provide approx- 



dw/dz as a function of redshift 
for 02,0 = 16, 10, 6, and 4 
(top to bottom) 
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Figure 5. Plot of dw/dz as a function of redshift. The values 
of 02,0 are listed in the order that the curves go top to bottom 
(02,0 = 16 is the red curve, 02, o = 10 is the orange curve, etc.), 
and chosen to be the same as those in Figure [3] 



imately 60 e-folds of inflation. The other values of 02, o were 
chosen to illustrate larger values of Swo- 

Figure [5] shows dw/dz, which, since w — Sw — 1, is the 
same as dSw/dz. This plot is consistent with what we might 
expect from Figure [3] it is only for z < 3 that dw/dz is 
significantly non-zero. 

We note here that others who have studied a variety of 
quintessence models near either maxima or minima in the 
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potential have obtained approximate analytical solutions for 
the evolution of w(a) (Dutta & Scherrer 2008, Dutta et al. 
2009, Chiba et al. 2009, Huang, Bond, & Kofman 2011). 
They have also concluded, as we will, that in these models 
w(a) is not well- fit by any linear function, including the 
popular toy model w(a) = wo+w a (l—a). Chiba in particular 
shows that the results of Dutta and Scherrer should apply 
to a quadratic potential such as the one we use. 



10 CONCLUSION 

In this paper, we have argued that the accelerated expansion 
in the early Universe (inflation) and the accelerated expan- 
sion today (dark energy) are caused by the same physical 
mechanism. If inflation is slow-roll, then dark energy should 
also be slow-roll, and this allows us a chance to observe a 
detectable difference 8w between w and —1. We have sug- 
gested that the simplest model of inflation, Linde's chaotic 
inflation, where the potential is quadratic, could also effec- 
tively model dark energy. For those who consider the super- 
Planckian values of the scalar field in this model problem- 
atic, we have suggested that the quadratic potential for in- 
flation and the second quadratic potential we use for dark 
energy could both be realized as N-flation describing many 
sub-Plan ckian fields that are closely clustered in mass. 

The idea that dark energy and inflation are caused 
by the same mechanism produces a real gain in predictive 
power. Since we know the minimum value of the scalar field 
for inflation to give the 60 e-folds of inflation we see in the 
observable Universe, we can, on the reasonable assumption 
that the two scalar fields would have similar initial values 
following a Gaussian distribution, estimate upper limit prob- 
abilities for observing different values of Swo = wq + 1 today. 

We also are able to use this slow-roll dark energy idea to 
make another prediction: that 8w oc 1/H 2 . Figure Ml shows, 
our numerical results seem to follow this formula. Our ana- 
lytical work shows that this is generically true in any slow- 
roll model of dark energy. This relationship can, we hope, 
provide a valuable template to observers seeking to look for 
a signal in noisy data. Indeed, if observers find this shape for 
w(z), it is a signature of slow-roll dark energy, which would 
be persuasive evidence that dark energy may be driven by 
a mechanism similar to that in inflation. It may also mo- 
tivate a focus on lower redshift observations, because this 
relationship implies that 5w will be close to zero for z > 3. 

Finally, a third result we derive is an approximate for- 
mula for the value of 8w today as a function of the value 
of the dark energy scalar field today, 02, o (equation (21 1). 
This formula is useful because if a deviation from w = — 1 is 



measured, it can be plugged in to tell us what the value of 
</>2,o is today in our model. This, in turn, may tell us some- 
thing about the initial value of the inflationary scalar field, 
as these fields are likely to have had initial values on the 
same order. Having an idea of the inflationary scalar field's 
initial value may be interesting because it can tell us how 
many e-folds more than the minimum 60 we observe the 
Universe as a whole underwent during inflation. 

Our model provides a simple physical picture of what is 
causing dark energy, and complements it with a simple ap- 
proximate formula for 5w = w + 1. We hope it will be useful 
as observers seek to interpret the upcoming data from BOSS 



(Baryon Oscillation Spectroscopic Survey) and Planck as 
well as other future dark energy missions. 
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